Diffusion and Transport Coefficients in Synthetic Opals 
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Opals are structures composed of the closed packing of spheres in the size range of nano-to-micro 
meter. They are sintered to create small necks at the points of contact. We have solved the 
diffusion problem in such structures. The relation between the diffusion coefficient and the termal 
and electrical conductivity makes possible to estimate the transport coefficients of opal structures. 
We estimate this changes as function of the neck size and the mean-free path of the carriers. The 
theory presented is also applicable to the diffusion problem in other periodic structures. 
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I. NANOSTRUCTURED THERMOELECTRICS 

A good thermoelectric material has low thermal con- 
ductivity K, high electrical conductivity a, and high See- 
beck coefficient, in order to maximize the thermoelectric 
figure of merit 



Z 



K 



(1) 



Z has units of inverse temperature and is generally 
quoted as ZT, with T the absolute temperature. For 
more than 40 years, the search for better thermoelectrics 
has not provided a material with ZT significantly larger 
than one. ZT of about 4 would make thermoelectric cool- 
ers capable to compete against the currently established 
devices based on gas-compression technology. 

A large effort has been done so far in or|der to improve 
the efficiency of thermoelectric materials Jd^El The search 
for better thermoelectrics can be classified by its length 
scale. The improvement was attempted in the micro- 
scopic, mesoscopic and macroscopic scale. In the micro- 
scopic scale, most of the search has been done by study- 
ing variations of the chemical composition of prospect 
materials and alloys. On the opposite range, the macro- 
scopic range involves the design of the thermoelectric de- 
vices. So far, most of the effort has been done in these 
length scales. Only recently, with the improvement of dif- 
ferent nano- fabrication techniques, the search for better 
thermoelectrics has been done in the mesoscopic length 
scale. This includes the study of superlattices, quantum 
dots and opals. It is interesting to mention an essen- 
tial difference between this and the previously mentioned 
length scales. The theoretical tools needed to evaluate 
the improvements in the micro and macro length scales 
are fairly well established. Both, the transport properties 
of compounds and alloys, and the macroscopic equations 
for device modeling are well known. However, this is 
not the case in the mesoscopic scale. In many possible 



structures, we are still working on the basic theoretical 
framework to understand and plan the new designs. A 
clear example of this is the discussion -around the ther- 
moelectric application of superlattices .a The purpose of 
the present work is to develope theoretical tools to eval- 
uate the transport coefficients of one of these mesoscopic 
structures: the synthetic opals. 

Synthetic opals are nanostructured materials com- 
posed of closed-packed spheres with a uniform radius a 
which can range from nanometers to microns. They are 
often made with glass spheres, but various techniques 
can be used to replace the glass spheres with almost any 
solid state material. The opals are usually sintered so 
that the points of contact between neighboring spheres 
become small necks. Currents of electricity or heat can 
go from sphere to sphere through these necks. Here we 
calculate the electrical and thermal conductivity of such 
opal structures. We ask the question: how does the opal 
structure reduce the electrical and thermal conductivity? 
Assuming that the Seebeck coefficient will not be affected 
by the opal structure, if the thermal conductivity is re- 
duced much more than the electrical conductivity, the 
opals could be useful thermoelectric materials. 

Currently, there is an experimental group at AUiedSig- 
nal Corporation preparing and developing these struc- 
tures to evaluate their potential as thermoelectrics. Our 
theory provides a theoretical tool to guide this evaluation 
and design.Ql3 The thermal conductivity of nanocompos- 
ites with regular structure have been studied by B. N. 
Bogomolov and collaborators arijd is, reported is a series 
of publication from 1995 to date.U til The object of these 
studies has been a Si02 nanocomposite opal with first 
order voids completely filled with NaCl. Their results 
are analyzed with a theory for composite materials de- 
veloped by R. E. Meredith and C. W. Tohiasli3 on the 
basis of the original work of Lord RayleighO. The the- 
ory is based in the solution of the Laplace equation for a 
regular array of spheres embedded in a different medium. 



This theory can not be appUed to the case for which the 
spheres are touching each other, which is the focus of our 
present interest. 

In Sec. m we relate the diffusion coefficient in the opal 
with the transport coefficients that are relevant for a ther- 
moelectric analysis. The diffusion coefficient is evaluated 
with a random- walk model as explained in Sec. |lll| and 
this model is analyzed and solved in Sees. [V and[V|. The 
results of applyin g th is model to an fee opal structure are 



discussed in Sec. VI 



II. 



TRANSPORT COEFFICIENTS AND 
DIFFUSION 



Given a bulk material with thermal conductivity k and 
electrical conductivity a, we want to estimate the ther- 
mal and electrical conductivities, k(°p) and (T*^°p\ of a 
synthetic opal made of spheres of this material. The ori- 
gin and magnitude of the reduction in the transport co- 
efficients depends on the mean free path of the carriers ^ 
compared with the diameter of the spheres d. li i « d 
the effect is due to a change in the boudary conditions, 
i.e. the thermal conductivity of the overall structure can 
be obtained by solving Laplace equation in the opal ge- 
ometry, li (. p^ d the boundary scattering produced by 
the sphere surface and the necks becomes important. In 
addition to the geometrical effect, which is also present 
in this case, there is a reduction of the intrinsic conduc- 
tivities due to the presence of the boundaries. In order 
to cover both situations with our model calculation, we 
will relate the transport coefficients to the diffusion coef- 
ficient of the opal and calculate this diffusion coefficient 
as a function of the mean free path and the opal geomet- 
ric parameters, which are the diameter of the spheres d 
and the distance between its centers a. 

The lattice thermal conductivity is proportional to the 
diffusion coefficient for phonons 



= CL», 



ph 



(2) 



where C is the specific heat of the material. Since we 
are mostly interested in semiconductors at room teipoer- 
ature, it is accurate to use the Einstein's relationlifl to 
relate the mobility to the diffusion coefficient for electri- 
cal carriers 



M 



keT 



(3) 



where e is the electron charge and ks is the Boltzmann 
constant. 

With Eqs. (|2[) and (H) in mind, to estimate the the- 
moelectric figure of merit of opals we have to solve the 
diffusion problem in the opal structure. In this way we 
will obtain a relation between the diffusion constant D 
of the bulk material and the overall diffusion constant of 
the opal D^°^\ The thermoelectric figure of merit Z is 
proportional to Dei/Dph- To evaluate the improvement 



produced by the opal we have to compare this ratio in 
bulk and in the opal. 



III. DIFFUSION COEFFICIENT AND RANDOM 
WALKS 

Consider a particle moving with velocity v between col- 
lisions, each collision randomly changing the direction of 
the movement. In this situation the diffusion coefficient 
is given by 






(4) 



where {r(t)'^) is the mean-square displacement of the par- 
ticle and t the total time of the random walk. This can 
be calculated by running a random walker that performs 
a step of lensth £ at every time step of length r. S. 
Chandrasekarlla solved the diffusion problem of a random 
walker in an homogeneous media and obtained 



(r2) = fN 



(5) 



where N is the number of steps. By replacing this ex- 
pression for (r^) in Eq. (|j) for the diffusion coefficient, 
we obtain 



D = -vl 
6 



(6) 



the 



where v is the velocity of the random walker, i.e 
ratio £/t. 

The expression given in Eq. (m for the mean-square 
displacement is only valid in bulk. If the carrier moves 
inside the spheres of the opal structure we have to cal- 
culate (r^) in a different way in order to use Eq. (0) and 
obtain the diffusion coefficient of the opal D'^°p\ 



IV. DIFFUSION COEFFICIENT IN OPALS 

We consider the carriers as classical particles perform- 
ing a diffusive motion inside the material that forms the 
opal. The particle will move inside a sphere crossing from 
time to time the neck to a neighboring sphere. If the 
length of the jumps is smaller than the diameter of the 
balls, the diffusion coefficient inside a ball will be that of 
the bulk material with a microstructure corresponding to 
grains of the size of the spheres. The effective diffusion 
coefficient of the opal structure will be determined by the 
diffusion from sphere to sphere. 

One method to evaluate this diffusion coefficient is the 
numerical simulation of the walker in the opal structure. 
However, if the diameter of the necks is small, as it is 
the case for the experimental situation, the walker will 
spend most of the time wandering inside a sphere and 



only sporadically crossing to the next sphere. As a con- 
sequence, most of the computer time will be lost in a dif- 
fusive movement inside the spheres and the diffusion co- 
efficient calculation, which is essentially related with this 
motion from sphere to sphere, will become extremely ex- 
pensive in computer time. A different method is needed. 

Here we present such a method. The central idea is to 
integrate out the motion of the carrier inside the sphere 
while keeping attention on the motion from sphere to 
sphere. The diffusion in opals, with this picture in mind, 
is similar to the diffusion of a particle in a lattice, the 
lattice made up by the spheres. The standard theory of 
diffusion in a lattice is due to Chudley and Elliot Jl3 How- 
ever, their theory was for a particle diffusing from site to 
site whithout spending any time in the site. The diffusion 
in an opal is different and the theory must include the 
feature that it takes time to diffuse within the spheres, 
from neck-to-neck. The opals are on a lattice, but the 
theory must include the diffusion within the spheres. 

Define 7ij(t) as the probability per unit time that a 
particle departs through the neck i at time t if it entered 
the sphere through neck j at i = 0. Define /i{Rn,,i) as 
the flux of particles departing a sphere centered at R„ 
through the neck labeled i at time t. Note that this quan- 
tity is not the net flux, which is the number of particles 
leaving minus the number entering, ft is just the num- 
ber leaving. The number entering the sphere through 
that neck is included as the flux leaving the neighboring 
sphere at R„ -1-5^ through neck i. Here we define the con- 
jugate to neck i as i. It is the label of the corresponding 
neck on the neighboring sphere. A neck i on one sphere 
is connected to i on the neighbor. 

The rate equations for particle motion are 

Jo 

+ S,sSnS{t) (7) 

The first term on the right is the outgoing flux due to 
all of the incoming particles at an earlier time. The last 
term on the right is the source term which starts the 
particle diffusion. The time integral can be eliminated 
by a Laplace transform 

/•oo 

E,{Rn,p)= / dte-P'MRn,t), 
Jo 

/"OO 

r„(p)=/ die"P*7^,(i). (8) 

After transforming Eq. (0) that describes the particle mo- 
tion, it becomes 

F,{Rn,p) = J2 ^^3{p)F]i'^n + 8j,p) + 5n (9) 



The lattice properties are taken into account by taking a 
Fourier transform 



G,(k,p)=^e^''-^"F,(R„,p), 



(10) 



which leads to the final form of the equation of motion 



G,(k,p) = ^r,-G,(k,p) 



,— ik-^. 



fll) 



where we have used 8j — —Sj. 



Eq. (nil) is a matrix equation which is solved for the 
functions Gi(k,p). However, the diffusion properties are 
determined by the properties of the matrix 



Mij = Sij - Tij cxp(-ik • Sj] 



(12) 



The poles of Gi are determined by the zeros of the de- 
terminant of this matrix. In taking the inverse Laplace 
transform, the time dependence is given by the poles of 
Gi, which again are the zeros of det|My|. The small- 
est pole is defined as that with the smallest value of p. 
The long time behavior (diffusion) is determined by the 
smallest pole. We now discuss some properties of the 
matrix M, and the propagators 7 that will be helpful 
in obtaining a form for the diffusion coefficient in opals 

£)(op). 

(i) If there is no absorption, then all of the particles which 
enter the sphere eventually must leave it sometime. 



1 = 



/•OO 



(0) 



(13) 



Symmetry shows that this result is independent of i. If 
we define 



TT{p) = J2r,,{p), 



the identity in Eq. (|13D, takes the form 

rT(o) = 1 . 



(14) 



(15) 



(li) // k == andp^O, then det \M\ =0. If k = each 
row of the matrix has the same elements, but in different 
order. As a consequence, the eigenvector ipo which has 
all elements equal to one has the feature that 



M4,o = il-TTipWo 



(16) 



Therefore at p = then MifjQ — 0. The vector ^0 has 
an eigenvalue of zero at p = 0. Since the determinant is 
the product of the eigenvalues, this also means that the 
determinant of Mij is zero in this case. 
(Hi) Atk^ det |M| has a factor 1 - Trip)- This is a 
direct consequence of the previous property. 

All of the above results are at k = 0. For diffusion one 
wants large values of R„ which means small values of k. 
The expansion of the determinant at small p and small 
k — |k| has the general form 



det \M\ = pu{p)U{p) + {kafV{p) 



(17) 



where U{p) and V{p) are polynomials in the propagators 
Tij(p), a is the nearest-neighbor distance, and u{p) = 
(1 — rT(p))/p- As a result, the diffusion coefficient of the 
opal is given by 



£)(.°P) — 



F(0) 



u(0) t/(0) ' 



(18) 



where m(0) is the limit of u{p) for p —^ 0. 

These general results are illustrated for some relevant 
common lattices. For each particular lattice, the propa- 
gator Tij depends only on the angle between necks i and 
j and we will use the notation Fg, where 9 is the angle 
between the normals to the neck surfaces. 



Simple Cubic Lattice 

The simple cubic lattice has four necks at an angle of 
90 degrees, as a consequence 



F, 



(sc). 



^{sc). 



(sc) 



T (p)=Fr^(p)+4Fgo 
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180 



(P) 



and the diffusion coefficient is given by 



D(op) 



180 



2F 



[sq) 
90 



6w(0) 2T^„^^ 



^(sq) 



Face Centered Cubic Lattice 



(23) 



(24) 



Linear Chain 

We first present the simple case of the linear chain in 
order to illustrate the method. In this case, each ball has 

two necks. We will define Fq (p) — Fii(p) = F22(p) as 
the Laplace transform of the probability that the walker 
enters and leaves through the same neck, and V^^q {p) = 
^i2{p) = ^2i{p) corresponding to the probability of the 
walkers leaving the sphere through the opposite side. In 
this case the matrix M has the form 



^('c) 



1 



^ 180^ 
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^^ •' c~ika 
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(19) 



The diffusion coefficient for this one dimensional opal will 
be 



Jjiop) 



F 
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180 



2u(o) r^''^) 



(20) 



Square Lattice 

In the case of the two dimensional square lattice we 
have three different possibilities, the walker can go back 
through the same entry neck, can go out through one 
of the two necks that form an angle of 90 degrees with 
the incoming trajectory, or can leave the sphere through 
the opposite neck. The first possibility is represented by 



^(s?)/ 



^(sq)i 



r};'^'ip), the second by F^7^(p) and the last by Vg^ip). 



^(sQ)! 



In this case 

risq). 



-^(sq) 



(sg), 



-^(sg)/ 



{p)^Tr>{p) + 2T^^,^>ip)+T\f>ip) 
and the diffusion coefficient is given by 



£){op) _ 



^{sq) 
180 



Asg) 
■90 



MO) rl,f + To 



(sg) 



(21) 



(22) 



The fee lattice corresponds to the close packing of 
spheres observed experimentally in opals. In the fee lat- 
tice each sphere is connected through twelve necks to 
its nearest neighbors. Assume that we enter the sphere 
through neck 1, in the fee lattice there are four necks 
with a normal forming an angle of 60° with the normal 
of neck 1, two at 90°, four at 120°, and one neck in the 
opposite side of neck 1 at 180°. 

The result for the fee lattice is: 



^ifcc). 



Afcc 



(fee), 



.(fee) 



[p) = T^r'ip) + ^^'^'ip) + ^nr'ip) + 



4F 



(fee) 
120 



-60 



(25) 



and 



]j(op) 



3u(0) 

yifcc) 
^ 60 



(26) 
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{o) + ni^^'{o) + sr\i^o>io) + r\Z'io) 



.(fee). 



.(fee) 



(fee). 



.(fee). 



(0)+3F^r(0) + F^r^(0) + Fi-^^(0) 

The final stage of our calculation of the diffusion coef- 
ficient in opals is the determination of the particle prop- 
agators inside the sphere, Jij{t), and this is the focus of 
the next section. 



V. PROPAGATORS INSIDE THE SPHERES 

The propagator inside the sphere, ^ijit), was defined 
as the probability of a particle leaving through neck j 
at time t if entered the sphere through neck i at t = 0. 
We obtain this propagators by numerical simulation of 
the walker. A specific behavior of the walker when it 
hits the surface has to be specified to do the simulation. 
We found that the final result is not strongly dependent 
on the particular scattering at the surface as long as the 
particles are not absorbed. All the results shown here 
have been performed assuming that the walker starts a 
new random step after hitting the wall. 

Since it is experimentally observed that the synthetic 
opal adopts a closed-packed arrangement of the spheres 



we will concentrate in the simulation of an fee arrange- 
ment of spheres. In Fig. |l| we show the time dependence 
of the propagators for an opal in which the diameter of 
the spheres d — 1.05 a and the step length of the walker 
£ — 0.06 a. The propagator for leaving through the same 
neck of entry, 70 (i), is the only one that does not go to 
zero for t —^ 0. There are many events in which the 
walker spends just a few steps inside the ball and goes 
back through the same neck. To leave through other 
neck the walker needs at least a number of steps equal 
to the ratio between the shortest distance between these 
necks and the step length. Of course this contribution 
corresponds to events with very low probability in which 
the walker reaches the other neck in an almost ballistic 
way. All the propagators converge to the same exponen- 
tial decay for long times. It is clear that after a long 
time of wandering inside the sphere, the walker looses all 
memory of the entry neck. 
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FIG. 1. Propagators as a function of time of a random 
walker inside an opal structure with diameter of the spheres 
d 5% larger than the intercenter spacing a. The mean free 
path £ = 0.06a. 

To calculate the diffusion coefficient we need the 
Laplace transform of this functions for p — that is the 
integral over all times, and the function u also evaluated 
a.t p — 0. From its definition we see that 



uiO) = - 



dTTJp) 
dp 



p=0 



t-iT{t)dt . 



(27) 



Both integrals are very easy to obtain from the numerical 
sinrulation. 



VI. DIFFUSION COEFFICIENT OF A FACE 
CENTERED CUBIC OPAL 

We have used the theory described above to evaluate 
the diffusion coefhcient of an fee opal structure. Figure 
shows the result as a function of the step length, for 
different sphere diameters. From the figure is clear the 
linear behavior of the diffusion coefficient with the step 



length, at least when it is shorter than 15% the sphere 
diameter. The diffusion coefficient of the bulk material 
is show as a reference. 
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FIG. 2. Diffusion coefRcient as a function of the mean 
free path for different diameters of the spheres. The diffusion 
coefficient is expressed in units of the velocity of the carriers 
times the lattice constant. The mean free path used for the 
plot is an average of the time step. The diffusion coefficient 
of the bulk material D — vl/Q is show as a reference. 

The result that the diffusion coefficient of the opal is 
a lineal function of the mean free path of the carriers 
is not promising for the opals as thermoelectric materi- 
als. This linearity implies that the ratio D^°i ID°^ will 
be approximately the same as the ratio Dei/Dph and no 
increase of the thermoelectric figure of merit is to be ex- 
pected. 

The slope depends on the diameter of the balls and 
defines an effective speed of the carrier. This dependence 
is shown in Fig. 0. This figure shows clearly that the 
diffusion coefficient goes to zero if the diameter of the 
spheres is smaller than the distance between the centers. 
This is because we assume a zero conductivity for the 
intersphere material. When the diameter is increased the 
effective velocity increases approaching that of the bulk 
material. 




d/a 

FIG. 3. Effective velocity as a function of the diameter of 
the spheres in the opal. 



VII. SUMMARY AND CONCLUSIONS 

We have developed a theory for diffusion of particles 
in an opal structure. With this theory, we calculate the 
diffusion coefficient of the opal as a function of the dif- 
fusion coefficient of the constituent material, the size of 
the spheres and the necks between them. We apply this 
theory to some of the most common lattices and provide 
an explicit expression for the diffusion coefficient in each 
case. 

The experimentally most relevant lattice is the fee. 
The application of this theory to a synthetic fee opal 
structure shows that the diffusion coefficient is linear with 
the mean free path of the carriers. This result holds up 
to mean free paths of the order of 15% the diameter of 
the spheres. The effect of the opal structure is to re- 
duce the effective speed of the carriers. We calculate the 
dependence of this effective velocity with the size of the 
spheres. 

The linear dependence of the diffusion coefficient with 
the mean free path indicates that the overall reduction 
for electrons and phonons will be the same. Assuming 
that the Seebeck coefficient is not affected by the opal 
structure, we do not expect any increase in the thermo- 
electric figure of merit of opals compared with the bulk 
material. A mayor effect of the opal structure is expected 
in the case of a mean free path for phonons larger or of 
the order of the sphere diameter. However, this will cor- 
respond to a material which is a good thermal conductor 
in bulk and the reduction should overcome this unfavor- 
able starting point. 

As a final comment, there is an interesting property of 
all the expressions obtained for the diffusion coefficient in 
the different lattices, Eqs. @, (|2|), (|2|), and @. In 
all these cases, adding the numerator and denominator 
rT(0) is obtained, i.e. the sum is equal to 1. Even though 
this relation is very suggestive we could not find a good 
explanation for it. 

The theory developed in this work is applicable, not 
only to synthetic opal structures but, can be used to solve 
the diffusion problem in any periodic structure in which 
the walker spends some time inside the unit cell. 
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